Macroscopic quantum electrodynamics and duality 
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We discuss under what conditions the duality between electric and magnetic fields is a valid 
symmetry of macroscopic quantum electrodynamics. It is shown that Maxwell's equations in the 
absence of free charges satisfy duality invariance on an operator level, whereas this is not true for 
Lorentz forces and atom-field couplings in general. We prove that derived quantities like Casimir 
forces, local-field corrected decay rates as well as van-der-Waals potentials are invariant with respect 
to a global exchange of electric and magnetic quantities. This exact symmetry can be used to deduce 
the physics of new configurations on the basis of already established ones. 
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In the past, studies of quantum electrodynamic (QED) 
phenomena have often been restricted to purely electric 
systems, because effects associated with magnetic prop- 
erties are considerably smaller for materials occurring in 
nature. Two developments have recently triggered an in- 
creased interest in such magnetic effects: The first was 
the suggestion [3| and subsequent fabrication Q of artifi- 
cial metamaterials with controllable electric permittivity 
e and magnetic permeability /i, where left-handed mate- 
rials (LHMs) with negative real parts of e and /i are of 
articular interest. As had already pointed out in 1968 
the basis vectors of an electromagnetic wave prop- 
agating inside such a medium form a left-handed triad, 
implying negative refraction. Motivated by the progress 
in metamaterial fabrication, researchers have intensively 
studied their potentials, leading to proposals of a perfect 
lens with sub-wavelength resolution Q as well as cloak- 
ing devices and predictions of an unusual behaviour of 
the decay of one or two atoms in the presence of LHMs 

Another, closely related motivation for considering 
magnetic systems was due to the fact that dispersion 
forces ^] have gained an increasing influence on micro- 
mechanical devices where they often lead to undesired 
effects such as stiction Q. The question naturally arose 
whether LHMs could be exploited to modify or even 
change the sign of dispersion forces. Forces on ex- 
cited systems might indeed be influenced by LHMs [lo| . 
Ground-state forces are not as easily manipulated be- 
cause they depend on the medium response at all fre- 
quencies, whereas the Kramers-Kronig relations imply 
that LHMs can only be realized in limited frequency 
windows. However, the controllable magnetic properties 
available in metamaterials can still have a large impact 
on dispersion forces: The dispersion forces between elec- 
tric and magnetic atoms [ll| or bodies 12 1 differ both 
in sign and power laws from those between only electric 
ones. On the search for repulsive dispersion forces, in- 
teractions of electric /magnetic atoms [l3], plates 3 isl ] 
and atoms with plates la, III have been studied; more 



complex problems such as atom-atom interactions in the 
presence of a magneto-electric bulk medium [31, plate 
[lit sphere [ISl have also been addressed. Reductions 
or even sign changes of the forces have been predicted 
for such scenarios and have been attributed primarily to 
large permeabilities rather than left-handed properties. 

Metamaterials have thus considerably increased the 
parameter space at one's disposal for manipulating QED 
phenomena. An efficient use of this new freedom re- 
quires the formulation of general statements of what 
might be achieved in principle. Working in this direction, 
upper bounds for the strength of attractive and repul- 
sive Casimir forces have been formulated [3] and it has 
been proven that the force between two mirr or-sy mmetric 
purely electric bodies is always attractive 2l|. In the 
present Letter, we establish another such general prin- 
ciple on the basis of the duality of Maxwell's e quat ions 
under an exchange of electric and magnetic fields 13, 3] , 
also known as electric/magnetic reciprocity within a gen- 
eralised framework of classical electrodynamics 2J]. In 
particle physics, duality has been discussed as a sym- 
metry of the M — 4 supersymmetric Yang-Mills theory 
(ist - We willprove its validity in the context of macro- 
scopic QED [all] and show that under certain conditions, 
quantities such as decay rates and dispersion forces are 
invariant with respect to a global exchange of electric and 
magnetic properties. The parameter space to be consid- 
ered in the search for optimal geometries and materials 
will thus be effectively halved. 

We begin by verifying duality for macroscopic QED 
in the absence of free charges and currents. We group 
the fields into dual pairs {y/enE, ^^H), {y/JIoD, ^Je^B) 
and {JJIqP, ^/sq^qM), so that Maxwell's equations read 
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Maxwell's equations are invariant under the general 
SO (2) duality transformation 



= v{e) 
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which may equivalently be expressed as a U(l) transfor- 
mation when introducing complex Riemann-Silberstein 
fields [2^. The invariance of Maxwell's equations under 
this rotation can be verified by multiplying Eqs. ([T|)-(l3]) 
by 'D{9) and using the fact that 'D{9) commutes with the 
symplectic matrix in Eq. Note that the grouping 

into dual pairs is solely due to the mathematical struc- 
ture of the equations and is in contrast to the fact that 
E, B and £), H are the pairs of physically corresponding 
quantities. 

For it to be a valid symmetry of the electromagnetic 
field, duality must also be consistent with the constitutive 
relations. In the presence of linear, local, isotropic, dis- 
persing and absorbing media, the constitutive relations 
in frequency space can be given as 
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where e — e{r,oj) and fj.~ij,{r,uj) denote the relative elec- 
tric permittivity and magnetic permeability of the media 
and J'n and A^n are the noise polarisation and magneti- 
sation which necessarily arise in the presence of absorp- 
tion. Invariance of the constitutive relations ^ under 
the duality transformation requires that 
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This condition is trivially fulfilled if e = /i (including 
both free space and the perfect lens, e — fj, — —l 0]), 
where duality is a continuous symmetry. For media with 
a non-trivial impedance, the condition ([6]) only holds for 
9 = n7r/2 with n € Z. The presence of such media thus 
reduces the continuous symmetry to a discrete symme- 
try with four distinct members, whose group structure 
is that of Z4. For 9 = n7r/2, Eqs. Q and © imply the 
transformations 

„2 a „;„2 , 
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E, H: 


E* 


= cjioH, 


H* = 


-E/{ciio) 


D, B. 


D* 


= ceo-B, 


B* = 


-n/iceo) 


P, M: 


P* 


= M/c, 


M* 


= cP 


Pa, Ma: 


Pa 


= Ma/c, 


MX = 


-- -cPa 


d, rh: 


d* 


— rh/c, 


rh* - 


- ~cd 


Pn, Mn: 


Pn - 


-- /iMn/c, 


MS = 


-cPn/e 


fe, fm'- 


n = - 


K^^/\^J■\)fm, 


fm ~~ 




e, n: 


e 


* = 




— e 


Q, /?: 


a* 


= /3/c^ 


13* 





TABLE I; Effect of the duality transformation with — -k/2. 

Maxwell's equations ((T|) and together with the 
constitutive relations (O for the electromagnetic field in 
the absence of free charges and currents, are thus invari- 
ant under the discrete duality transformations = n7r/2, 
n e Z given by Eqs. (g]), ^ and jS]). This is not only 
true for the equations of motion, but clearly must also 
hold on a Hamiltonian level. To see this explicitly, re- 
call that the Hamiltonian of the medium-assisted field is 
given by = Ea^g,™ / d^'' /o°° dcj /|(r, tj) •/A(r, 

where the fundamental bosonic fields f\ are related 
to the noise terms via 
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Combining Eqs. JT]), ^ and 
damental fields transform as 



(9) 

one finds that the fun- 
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for 9 — mi/2, so that Hp — Hp. It is sufficient to focus on 
the single duality transformation 9 = tt/2 as summarised 
m Tab. H which is a generator of the whole group. 

Let us next turn our attention to Lorentz forces and 
the coupling of the medium-assisted field to charged par- 
ticles: We recall that the operator Lorentz force on a 
neutral body occupying a volume V can be given as 



dA-{eoE{r)E{r) + —B{r)B{r) 
dv I /^o 



eoE^r) + —B\r) 



~ea-J^d^rE{r)xB{r) (11) 

(I: unit tensor), while that on a neutral atom with po- 
larisation Pa and magnetisation Ma reads [^, 27 1 



F = Va I dV 



PAir)xrA-B{r) 



PA{r)-E{r) + MA{r)-B{r 
d 



At 



d-^rPA{r)xB{r). (12) 
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The coupling of one or more atoms to the medium- 
assisted electromagnetic field can in the multipolar cou- 
pling scheme be implemented via [1, 



Haf = - 



PA(r)-E{r) + MA{r)-B{r) 

+ m-/PA{r)xpA-B{r)\, (13) 



when neglecting diamagnetic interactions. Using the 
transformation behaviour given in Tab. [H it is imme- 
diately clear that neither the Lorentz forces on bodies 
or atoms nor the atom-field interactions are duality in- 
variant on an operator level. Even for atoms and bodies 
at rest with time-independent fields, duality invariance 
is prohibited by the unavoidable noise polarisation and 
magnetisation in the constitutive relations (O. 

That said, we will show that effective quantities de- 
rived from the above operator Lorentz forces and atom- 
field couplings do obey duality invariance when consid- 
ering atoms and bodies at rest and not embedded in a 
medium. In particular, we will consider the Casimir force 

M 

F^-- r Aij dA.{G(i)(r,r,iO + G(l(r,r,iO 
- \ Tr [g« (r, r, i^ + G«„(r, r, i^] l} , (14) 
the single- and two-atom vdW potentials 0, [13, Izst 
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(a, /3: atomic polarisability, magnetisability) and the 
atomic decay rate |29[ 
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(|n): initial atomic state, ujnk' atomic transition fre- 
quencies; dnk, mnk- electric, magnetic dipole matrix el- 
ements). Here, G*-^-* is the scattering part of the classical 
Green tensor, where a left index e, m indicates that G 
is multiplied by \Lojc = — C/c or Vx from the left and 
a right index e, m denotes multiplication with iiojc^ 

—£,/c or xV' from the right. The Casimir force and the 
single-atom vdW force are the ground-state averages of 
the above operator Lorentz forces, while the atomic po- 
tentials and rates follow from the atom-field coupling. 

To prove the duality invariance of the above quantities 
(fT4|) -(fT7 |) . we note that the Casimir force depends solely 
on the classical Green tensor 



Vx 



■ Vx -~^e{r,uj) 



G(r,r',a;) = 5{i — r' 



(18) 

while vdW forces and decay rates also depend on a, (3, 
d and rh. While the transformation behaviour of the 
latter quantities under duality follows immediately from 
that of e, /i. Pa and Ma (see Tab. the transformed 
Green tensor, which is the solution to Eq. (|18p with e 
and fjL exchanged, can be determined as follows: We first 
note that Maxwell's equations ([1]), ([2]) together with the 
constitutive relations ([5]) are uniquely solved by [ij 

E{r,u:) = -So' I d\'Gee{r,r',u;)-P^{r',u;) 
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d^r' Gemir,r',uj)-Mj^{r\uj) 
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B{r,uj) = ~cfio J d^r'G,neir,r',uj)-P^{r',Lj) 
-Mo / d\'G,nmir,r',Lu)-M.Nir',i^), (20) 



D(r,uj) -- 
dV 



eir,uj 



dh-' Gera{r,r',uj)-Mj^{r',uj) 



e{r,uj)Gee{r,r',uj) - 6{r-r') 



(21) 



d^r' Gme{r,r\uj)-Pj^{r',uj) 



dV 



Grarn{r,r',Uj) 



+ S{r-r') 



■M^{r',u). (22) 



The invariance of Maxwell's equations implies that this 
solution remains valid after applying the duality trans- 
formation. Taking duality transforms of Eqs. (jl9p and 
pop , the unknown transformed Green tensor appears on 
the rhs of these equations, whereas the transformations 
of all other quantities occurring in the equations can be 
determined with the aid of Tab. [J After using Eqs. (fTO|) - 
to express the resulting fields on the Ihs in terms of 
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s-iid and equating coefficients, one obtains the 
following transformation rules: 

G*e(r,r',w) = ^i'^{r,uj)Gmni{r,r',uj)n-^{r',uj) 

+ fi~\r,uj)6{r-r'), (23) 

Gem(^>'f"'.^) = - ^i^\r,uj)G,ne{r,r',uj)e{r',uj), (24) 

G;;,,(r,r',w) =-£(r,c^)Ge™(r,r',c^)Ai-i(r',w), (25) 
^limir,r',uj) = e{r,uj)Geeir,r',uj)e{r\uj) 

- e{r,Lu)S{r-r'). (26) 

The duality invariance of dispersion forces and decay 
rates follows immediately. Using Eqs. (^5]) and and 
noting that the S function does not contribute to the 
scattering part of the Green tensor, it is seen that the 
Casimir force on a body is unchanged when globally 
exchanging e and /i, provided that the body is located 
in free space. The duality invariance of the vdW poten- 
tials p5|) and (fTH)) also follows from the transformation 
rules This invariance with respect to a simul- 

taneous exchange e ^ /i and a again only holds 

if s(rj^/B) = /^('Ta/b) = 1- In contrast to the Casimir force, 
this does not mean that the atom has to be located in 
vacuum, but merely implies that for atoms embedded in 
media, local-field corrections must be included via the 
real-cavity model in order to insure invariance [30| . 

Duality invariance can be used to obtain the full func- 
tional dependence of dispersion forces in given scenar- 
ios on the atomic and medium parameters from knowl- 
edge of the respective dual scenario. For instance, it 
has recently been shown that in the retarded limit the 
vdW potential of two polarisable atoms reads U{rAB) = 

''^As] when includ- 
ing local-field corrections [30|]. Making the replacements 
a — > £ ^ /i, one can immediately infer UItab) — 

^ab] magneti- 

sable atoms. The utility of this principle becomes even 
more apparent for complex problems like the interaction 
of two atoms in the presence of a magneto-electric object 



3, 



20|. Finally, using the fact that two purely electric, 
mirror-symmetric bodies always attract |2l| , we can im- 
mediately conclude that so do two purely magnetic ones. 

In addition, Eqs. ((23|) and ((26|) imply the duality in- 
variance of the decay rate p?)) since the i5 functions do 
not contribute to the imaginary part of the Green ten- 
sor; again, local-field corrections have to be included for 
atoms embedded in media. This symmetry can be ex- 
ploited, e.g., to obtain magnetically driven spin-flip rates 
of atoms in specific environments from known electric- 
dipole driven decay rates. 

In conclusion, we have shown that dispersion forces on 
atoms and bodies as well as decay rates of atoms are 



duality invariant, provided that the bodies are located 
in free space at rest and that local-field corrections are 
taken into account when considering (stationary) atoms 
embedded in a medium. The established symmetry oper- 
ation of globally exchanging electric and magnetic body 
and atom properties is a powerful tool for obtaining new 
results on the basis of already established ones. The in- 
variance can easily be extended to other effective quanti- 
ties of macroscopic QED such as frequency shifts, heating 
rates or energy transfer rates. 
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